Interference effects in isolated Josephson junction arrays with geometric symmetries 



O 

o 

(N 
cn 



o 
o 

u 

Oh' 



X3 
O 

o 



> 

m 
(N 
(N 
O 



I 

O 

o 



X 
S3 



Dmitri A. Ivanov", Lev B. loffe^, Vadim B. Geshkenbein" , and Gianni Blatter" 
" Theoretische Physik, ETH-Honggerberg, CH-8093 Zurich, Switzerland 
Center for Materials Theory, Physics Department, Rutgers University, Piscataway, NJ 08854 USA 

(February 9, 2001) 

As the size of a Josephson junction is reduced, charging effects become important and the supercon- 
ducting phase across the link turns into a periodic quantum variable. Isolated Josephson junction 
arrays are described in terms of such periodic quantum variables and thus exhibit pronounced 
quantum interference effects arising from paths with different winding numbers (Aharonov-Casher 
effects). These interference effects have strong implications for the excitation spectrum of the array 
which are relevant in applications of superconducting junction arrays for quantum computing. The 
interference effects are most pronounced in arrays composed of identical junctions and possessing ge- 
ometric symmetries; they may be controlled by either external gate potentials or by adding/removing 
charge to/from the array. Here we consider a loop of N identical junctions encircling one half super- 
conducting quantum of magnetic flux. In this system, the ground state is found to be non-degenerate 
if the total number of Cooper pairs on the array is divisible by N, and doubly degenerate otherwise 
(after the stray charges are compensated by the gate voltages). 



I. INTRODUCTION 

Josephson junction arrays are excellent tools for ex- 
ploring quantuntmechanical behavior in a wide range of 
parameter spaceB. The charge and phase on each island 
provide a set of conjugate quantum variables, allowing for 
dual descriptions of the array either in terms of charges 
(Cooper pairs) hopping between the islands or in terms of 
vortices hopping between the plaquets of the arrayB. This 
opens the door for the set up and manipulation of inter- 
esting quantum interference effects in Josephson junction 
arrays: in a magnetic field, charges pick up additional 
Aharonov-Bohm phases and hence the properties of the 
array depend on the field strength. In the dual language, 
vortices moving around islands gain phases proportional 
to the gwierage charges on the islands (Aharonov-Casher 
phases )0^D. These features have recently been used in var- 
ious proposals for solid state implementations of qubits 
for quantum computing, based on either the chargeQ or 
phaseO □ degree of freedom in Josephson junction arrays. 

In this paper we study an interference effect in elec- 
trically isolated Josephson junction arrays which renders 
the ground state and excitation spectrum sensitive to the 
total charge on the array. This effect combines the dual 
descriptions in terms of charge or phase: i) in the limit 
where the charging energy Ec is much larger than the 
Josephson energy Ej, the fluctuations of charge on the 
islands are small and the system is equivalent to strongly 
repulsive bosons (Cooper pairs) hopping between islands; 
the total charge then determines the number of such 
bosons and hence the structure of the spectrum, ii) In the 
opposite limit of large Josephson energy Ej, the phase 
fluctuations are small and the spectrum is determined by 
phase tunneling between classical minima of the Joseph- 
son energy; the total charge on the array then determines 
the interference between different tunneling trajectories 



via the Aharonov-Casher effect. These interference ef- 
fects, where the ground state and the excitations may 
change degeneracies depending on the total charge on the 
array, are most pronounced in small arrays with geomet- 
ric symmetries; for the symmetric loops considered in this 
paper they coincide in the charge- and phase-dominated 
limits and persist at arbitrary ratio of the Josephson- to 
charging energy. In mathematical terms, the total charge 
on the array enters its symmetry group producing a cen- 
tral cxtcnsiorn. The irreducible representations are then 
classified by the total charge, and their dimensions (and 
associated level degeneracies) depend on the charge on 
the array. 

In experimental realizations of Josephson junction ar- 
rays, the charges and potentials on the islands are af- 
fected by a number of external factors such as differences 




FIG. 1. A superconducting phase qubit consists of a small 
inductance loop made from A'' superconducting islands con- 
nected by Josephson junctions and capacitively coupled to 
each other and to the ground (the case A = 3 is shown here) . 
The loop is placed in an external magnetic field producing 
one half superconducting flux quantum through the loop. 
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in island size, background charges, charges localized at 
impurities, etc. All these factors produce random charge 
offsets on the islands which distort the quantum inter- 
ference pattern (see e.g. the discussion of offset charges 
m Ref. 0). In order to observe the quantum interference 
effects predicted in this paper the potentials of the is- 
lands should be adjustable via external capacitive gates. 
The offset charges should then be compensated by ap- 
propriate tuning of the gate voltages; such a procedure 
was successfully carried out in a transport measurement 
involving two islandgli3 and should also be possible for 
the system considered in this paper. We also assume the 
absence of quasiparticles on the islands. Firstly, this re- 
quires the temperature to be sufficiently low to eliminate 
thermal quasiparticles. Secondly, the total number of 
electrons on the array must be even in order to exclude 
unpaired charge due to the parity effect. 

In the following we study the simplest setup exhibiting 
these interference effects, a closed low-inductance loop 
with TV > 3 identical Josephson junctions and pierced by 
one-half superconducting flux quantum (Fig. Such 
loops have been proposed as possible |j;ealizations for su- 
perconducting phase quantum bitsQ H and a successful 
quantum superposition of the two qubit states has been 
reported recentlyO. While in the original qubit design 
the loop has been made asymmetric in ordet,tp suppress 
the interference of different tunneling pathma, it is this 
interference that we study in the present paper where 
we assume the islands and contacts to be nearly iden- 
tical (the required precision will be specified at the end 
of the paper). The ground state degeneracy of such a 
symmetric Josephson-junction loop is found to depend 
periodically on the total charge Q with period 2eiV, N 
is the number of islands in the loop. In particular, in 
the phase-dominated regime Ej ^ Ec the Josephson 
energy of the fully frustrated loop exhibits two equiva- 
lent minima describing states with currents circulating 
in opposite directions. Tunneling between these minima 
produces a splitting A between these levels which is de- 
termined by the charges qi (in units of 2e) induced on 
the islands, 



j^p2lTi{qi+...+qN-l) ^ 



and which vanishes for (Q/2e) ^ (mod N) (the sum of 
all induced charges qi equals Q, and the expression (|l|) 
is symmetric under circular permutation of the islands). 
This result is the consequence of the interference of dif- 
ferent tunneling paths connecting the two minima where 
the relative phases of the tunneling amplitudes depend 
on the charges qi (a similar effect was predicted for the 
S-S-S double junction in Ref. |l^). The induced charges 
qi may be tuned by either gate voltages (redistributing 
the existing charge) or by adding/removing extra charge 
to/from the array. E.g. suppose that all charges qi have 
been set to zero by fine tuning the gate voltages, thus 



maximizing the splitting A. Then adding (or remov- 
ing) a charge Q at fixed relative potentials adds a value 
Q/2eN to each induced charge qi and closes the gap A: 
the splitting is present when the charge Qjle is divisi- 
ble by N and absent otherwise. A similar result applies 
for the opposite limit Ec S> Ej where the ground state 
is non-degenerate with the next excited state an energy 
Ec away if Q/2e is divisible by N ('insulating' state), 
while it is degenerate with an excitation gap of order Ej 
if Q/2e ^ (mod N) ('metalhc' state). Indeed, this 
Q-dependence of the ground state degeneracy will be ex- 
plained using symmetry considerations valid in the entire 
range of couplings Ej / Ec ■ 

The paper is organized as follows. In Section II, we 
define the model and prove the periodic dependence of 
the excitation spectrum on the total charge Q, for a sym- 
metric loop. In Sections III and IV, we treat the limits 
Ej ^ Ec and Ej <^ Ec, respectively. Section V con- 
tains the analysis of the symmetries of the loop. Finally, 
in Section VI we discuss the physical requirements for 
observing the charge-dependent interference effects. 



II. MODEL AND PERIODIC g-DEPENDENCE 
OF THE SPECTRUM 

The Hamiltonian describing the qubit in Fig. |l| takes 
the form 



-^Ui{ip^+i ~ ifii - ai^^+i). (2) 



Here, ipi are the phases on the islands, Qi = —id/dipi 
are the charge operators conjugate to (pi (we measure Q 
in units of 2e from now on), (C^^)y is the inverse ca- 
pacitance matrix, Vi are the gate voltages applied to the 
islands, Ui are the Josephson energies of the junctions, 
and ai_i_|-i is the electromagnetic vector potential induced 
by the external magnetic field (« -f 1 in the indices should 
be understood modulo N). This Hamiltonian acts on 
the wave functions ^'(vi, . . . , <fN) which are periodic in 
all their variables, 



■2tt,...,lpn) 
= ^'((^i,...,(^i,. 



(3) 



In general, the boundary conditions (^ may contain ar- 
bitrary phase shifts e*^' incorporating the effect of back- 
ground charges. They can be manipulated by the gate 
voltages Vi and we assume them to vanish through ap- 
propriate fine tuningl3 (this is equivalent to adjusting 
the zero positions of the gate voltages). We assume that 
this tuning is performed once in the beginning of the ex- 
periment and later measure the gate voltages relative to 
these reference values. 
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Since the potential term in the Hamiltonian (g) con- 
tains only phase differences Lpi — (fj, it has the sym- 
metry of rotating all the phases by the same angle: 
^(iy9i, . . . , (^jv) ^ '^{ifi + 5lp, . . . .i^N + Sip). Equiva- 
lently, the total charge Q — J2iQi conserved (i.e., 
the charge Q commutes with the Hamiltonian H) and 
we may project the Hilbert space onto the subspace with 
a given total charge Q before diagonalization, implying 
the transformation rule 

■f{ipi +S(p,...,(pN + S(p) = e'*'^'5*((^i, ...,(Pn) (4) 

for the simultaneous rotation of all the phases by Sip. 
From the periodicity of it immediately follows that 
Q is integer, i.e., the total charge must be a multiple of 
2e (this is in fact an implicit assumption when writing 
the Hamiltonian (^) in terms of phases only). In the fol- 
lowing we shall discuss the symmetric loop, postponing 
the effects of asymmetry till the end of the paper. Here, 
by symmetry we mean {C~^)ij = {C~^)i+k,j+k for the 
Coulomb term and equality of all Josephson terms Ui{ip). 
This implies that the loop is invariant under circular per- 
mutation of the islands. Also, if the flux through the loop 
is exactly one half flux quantum, the loop is invariant un- 
der 'flips' changing the direction of the current. 

The excitation spectrum of the symmetric loop peri- 
odically depends on the total charge Q with the period 
N (up to overall shifts): the unitary operator 

U -.^^ ^e'(¥'i+-+Viv) (5) 

increases the charge on all the islands by one, U~^QiU = 
Qi + 1, and therefore the sector with total charge Q maps 
onto the one with total charge Q + N. On the other hand, 

u-'Hu = H+ j2ic-'hQj + 1 E(^"')^^ + E 

i ij i 

If all the islands are equivalent, then the sum X]i(C'~^)y 
is independent of j and 

U-'HU = H+^{q + ^\+Y,V., (6) 

^ i 

where C = N[Y,,{C-^\jY^ = Y.^o'^^3 is the total 
capacitance of the loop. Thus the operator U maps 
the eigenfunctions of the Hamiltonian in the sector with 
charge Q onto eigenfunctions in the sector with charge 
Q + N shifting them in energy by a constant as given 
by (^j thus proving our statement about the periodic 
Q-dependence of the excitation spectrum. 



III. LEVEL SPLITTING IN THE E., > Ec LIMIT 

In the phase-dominated regime with Ej 3> Ec the 
low-energy states of the Josephson junction loop are de- 
termined by the classical minima of the Josephson en- 
ergy as corrected by weak quantum tunneling (due to 



(a) (b) 




FIG. 2. (a) One flux quantum moving around the j-th 
island; (b) One of the A'' possible tunneling trajectories: 
one flux quantum enters the loop through one of the junc- 
tions and leaves divided uniformly between all junctions; (c) 
The phase difference between two tunneling trajectories with 
flux entry through the junctions k\ and k2 is given by the 
Aharonov-Casher phase arising from one flux unit circling 
the (charged) islands in between. Note that closing the flux 
trajectory in (c) on the opposite side produces the same phase 
because Q is an integer. 

the finite charging energy). Here, we consider a sym- 
metric loop pierced by half a quantum of magnetic flux, 
'Ylii = ""i E^nd we choose to work in a gauge with 

<ii,i+i — tt/N for all i. In the limit Ej ^ Eq the only 
constraint on the potential Ui{ip) is the double degener- 
acy of the total Josephson energy as a function of the 
phases ipi (e.g., this requirement is satisfied for tunnel- 
ing junctions with Ui{ip) — —Ej cos ip and > 3). The 
two potential minima are determined by the phase con- 
figurations ipi = and ipi = {2Tr/N)i (and all configura- 
tions obtained from these two by the continuous symme- 
try ipi i—f ipi + Sip). These two minima involve different 
directions of the Josephson current, circulating the loop 
clockwise or counter-clockwise. The continuous symme- 
try ipi ipi + Sip, i = 1, . . . , N implies the quantization 
of the total charge Q as discussed above. 

In the following it will be convenient to incorporate the 
voltage terms ViQi into the quadratic kinetic term via a 
shift Qi Qi — Qi, where 

'-Eft^fv.-S^j+I (7) 

are the mean charges induced on the islands (we also 
include the i-independent contribution from the total 
charge Q for further convenience). This gauge trans- 
formation eliminates the ViQi terms in the Hamiltonian 
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(and adds a constant energy shift which we disregard). 
Under this transformation the wave functions pick up 
additional phases = exp(— q^i^J^'; the new 

wave functions ^' satisfy the twisted boundary conditions 

§((/?!,..., (Pj- +27r,...,(^A,) (8) 
= exp{-2Triqj)^{ipi, ...,(pj,..., ipn) 

and the invariance under the continuous symmetry Lpi i— > 
tfi + 6lp takes the form (c.f. (Q)), 

^{ipi +5lp,...,lpn + 5^)^ ^(<y5i, . . . , ^pn)- (9) 

We introduce the new variables 0+ ~ [ipi + . . . + (Pn)/N 
(conjugate to Q) and the gauge invariant phase drops 
over the junctions (f>i = ^Pi+i — 'Pi — <ii,i+i (out of N vari- 
ables only A^ — 1 are independent, since (f>i — — tt). 
The Hamiltonian then decouples into two, one involving 
only 0+ and the other only the phase drops (pi. The con- 
tinuous symmetry (^) implies that ^ is independent of 
(/>+, i.e., \1/ = ^'(01, . . . , 0jv)- The boundary conditions 
on ^ derived from (0) then take the form 



*('/'i, • • • , 0i + 27r, (pj+i - 27r, . . . , 

= exp(27rigj)*((?!)i, . . . , 0^, (/i^+i, . 



(10) 
) 



and can be given a simple physical interpretation (see 
Fig. H(a)): on moving a flux 27r (one flux quantum) 
around the j-th island the wave function gains a phase 
2Trqj . 

We are now prepared for the calculation of the tun- 
neling amplitude connecting the two semiclassical min- 
ima. Due to the equivalence of phase drops (pi dif- 
fering by a multiple of 27r, the coordinate space is a 
{N — l)-dimensional torus. After unfolding the torus 
onto the plane, each point on the torus is represented 
by a lattice of points on the plane. We choose the point 
{(pi = —tt/N} as the representative of one of the two 
minima. The N neighboring minima with coordinates 
{0, = n/N, i ^ k; (pk ^ n/N - 2tt}, fc = 1, . . . , TV then 
correspond to the other minimum. Accordingly, there 
are N different optimal tunneling trajectories connect- 
ing the two minima and we have to add up the vari- 
ous amplitudes coherently in order to find the hopping 
amplitude. Inspection of the particular phase changes 
along each of these trajectories shows that they may 
be thought of as describing a flux 27r entering the loop 
through the junction k and leaving the loop equally dis- 
tributed (i.e., with a fraction 2Tr/N of flux per junction) 
among all the junctions, see Fig. ||(b). These trajectories 
exhibit different phase windings and thus are topologi- 
cally different. The relative phases between these tra- 
jectories may then be read off the boundary conditions 
for the wave function ^E*; e.g., the phase difference be- 
tween the tunneling amplitudes for the trajectories ki 
and ^2 equals the phase difference of the wave function 
^ at the points {(pi = tt/A^, i ^ ki; (pki = t:/N ~ 27r} 



and {(pi = tt/N i ^ k2\ (t>k2 = tt/A^ — 27r}. This phase 
difference is 27r(gfci+i + qkx+2 + • • • + fe) and may be 
interpreted as the phase generated by the unit flux en- 
tering the loop through junction fci and leaving through 
junction /c2 (see Fig. ^(c)). From this argument imme- 
diately follows the result (|l|) for the level splitting (the 
level splitting in the double-well problem is proportional 
to the absolute value of the tunneling amplitude). 



IV. LEVEL SPLITTING IN THE E., < Ec LIMIT 

In the charge-dominated limit Ej <C E^-At is conve- 
nient to work in the charge representationll3, where the 
operators Qi are diagonal and the Josephson term in the 
Hamiltonian (|^) takes the form (we restrict the discussion 
to tunneling junctions with Ui{ip) = —Ej cos{lp)) 



Hi 



Ej 



E 



in /r, 



h.c. 



(11) 



with {L^) the charge raising (lowering) operators on 
the «-th island, L'^\Qi) = IQi ± 1). The periodicity in Q 
allows us to restrict our analysis to the N charge sectors 
Q = 0, . . . , N — 1. Below we first neglect the coupling 
Ej and find the ground states of the Coulomb part of 
the Hamiltonian; hopping between the islands is then 
perturbatively included in a second step. 

We start with a diagonal matrix (C^^)ij and ignore 
the capacitive coupling of the junctions, i.e., the islands 
are coupled only to the ground but not to each other. In 
that case, the ground state of the Coulomb part of the 
Hamihonian is C^-fold degenerate (C^ = m/Q\{N-Qy. 
enumerates the number of ways to distribute Q parti- 
cles among N sites without double occupancy). Second, 
we solve the Hamiltonian (|ll|) projected onto this C^- 
dimensional subspace in order to find the level splitting 
at finite Ej. This is easily done by observing the equiva- 
lence of this problem to the tight-binding model for hard- 
core bosons on the circle with N sites. Mapping to free 





FIG. 3. The tight-binding spectrum on the circle with pe- 
riodic (a) and anti-periodic (b) boundary conditions applying 
for even and odd Q, respectively. The solid circles denote 
filled states at the bottom of the band. For Q not divisible 
by A'', the ground state is doubly degenerate. Q divisible by 
N corresponds to an empty (or filled) band. 
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fermions and taking into account the flux through the 
loop and the boundary conditions on the circle, one finds 
that the projected Hj describes a tight binding model 
for Q free fermions on a circle with N sites and peri- 
odic boundary conditions for even Q, while anti-periodic 
boundary conditions apply for odd Q (Fig. |^). This im- 
plies that the ground state is non-degenerate with a gap 
of order Eq if Q is divisible by N (the 'insulating state' 
with an empty band) and doubly degenerate (with the 
lowest excitation at an energy of the order of the intra- 
band level spacing Ej/N above the degenerate ground 
state) otherwise. 

Turning on the junction capacitance makes the parti- 
cles (Cooper pairs) repel each other and they tend to ar- 
range in configurations with maximal separation between 
them. The number of such configurations is generally less 
(or equal) than C^. We conjecture that even in this case 
the ground state level is degenerate unless Q is divisible 
by N . We do not have a rigorous proof of this statement 
but have verified it for N = 3,4,5 (in fact, for = 3 
the off-diagonal elements of the matrix (C~^)ij may be 
incorporated in a constant term oc and do not change 
the properties of the system). The degeneracy may, of 
course, be explicitly verified for any given N and Q. 

V. SYMMETRY ANALYSIS 

In the previous sections, we have seen that the ground- 
state degeneracy depends on the divisibility of the total 
charge Q by the number of islands and coincides in 
the two limits Ej <C Ec and Ej ^ Eq- In this section 
we show that this degeneracy is a consequence of the 
geometric symmetries of the loop and remains exact be- 
yond the perturbation theory around these limiting cases. 
We shall classify the irreducible representations of the 
symmetry group of the loop, and identify the representa- 
tions corresponding to the ground state. In both limits, 
the ground states correspond to the same representation, 
which implies that the degeneracy is also preserved in 
the intermediate parameter range, unless level-crossing 
occurs. We verify the absence of level-crossing in the 




FIG. 4. Geometric symmetries of the Josephson junction 
loop for odd and even A'' involving translations T and reflec- 
tions R. For even A'^, the two classes of reflections are related 
by R* = TR. 



several simplest cases numerically by exact diagonaliza- 
tion. 

The level degeneracies for general values of Ej /Eq are 
determined by the symmetry group of the Josephson- 
junction loop. The geometric symmetries are described 
by the dihedral group Djsi consisting of cyclic permuta- 
tions T of the islands (preserving the loop orientation; 
these are equivalent to a iV-fold rotation axis) and of re- 
flections R about diameters of the loop drawn through 
islands or through junctions, see Fig. ^. The reflections R 
reverse the flux through the loop and belong to its sym- 
metry group for flux zero or half-quantum. The dihedral 
group Diq involves symmetry operations; it may be 
characterized with the defining relations 

= 1, (12) 
i?2 = 1, (13) 
{TKf = 1. (14) 

These relations are obeyed by the operators represent- 
ing R and T at zero magnetic fiux through the loop. 
On the other hand, for the case of half a flux quan- 
tum piercing the loop the symmetry operators preserving 
the Hamiltonian need to be supplemented by additional 
gauge transformations. As a consequence, the above re- 
lations are modified, with the appearance of additional 
phase shifts. Explicitly, the operators corresponding to 
rotations T and reflections R have the form 

T5'((pi, . . . , ipN) = *(¥'2, ■ • ■ , (^Af, •^l), (15) 
R^{ipi,. ■■,ipN) = ^('PAr, ^N~l + 2tt/N, (16) 

ipi+2Tr{N ~l)/N), 

where the shifts {2,11 /N)k of the phase tpN-k compensate 
for the external vector potential in the Hamiltonian (|^) 
and guarantee its invariance. This modifies the relation 
(|l3| ) producing an additional phase shift: 

i?2 = exp(-27riQ/Ar). (17) 

The simplest way to make use of symmetry arguments 
is for the case A^ even and Q odd: we show that all states 
are degenerate in this case. Let us assume the opposite, 
i.e., that the eigenstate \'^) is non-degenerate. Acting 
with T and i? on j^f) we reproduce the state up to phases t 
and r which satisfy the relations = {tr)^ = 1 and = 
exp{—2TTiQ/N); this set of equations is inconsistent for 
N even and Q odd, hence all levels indeed are degenerate. 

To determine the level degeneracy and to treat the case 
of arbitrary A'^ and Q, we classify the representations of 
the symmetry group. To take into account the phase 
shift in (|l7|), we include such phase shifts as new central 
elements {^"|^=i} in the group. The resulting set of 
defining relations is 

= = (Ti?)2 = 1, R^ = Z-\ 
ZR = RZ, ZT = TZ (18) 
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TABLE I. Character table of the dihedral group Dm for 
odd: the two one-dimensional representations are de- 
noted by D'*'; the (A'^ — l)/2 two-dimensional represen- 
tations D**'' involve the roots = exp{2nifi/N), with 
M = l,...,(A-l)/2. 

{E} [T.^T-.}iN-i)/2 {T-^RCzI} 

D^^'l 1 1 ±1 

D^"' I 2 f^+t;," 



TABLE IL Character table of ZDn for A odd. The rep- 
resentations are classified by the charge Q which relates to 
the values of the central element Z — exp(27riQ/A). The 2A 
one-dimensional representations are denoted by the 
A(A^ — l)/2 two-dimensional representations D'^'"'" involve 
the roots — exp{2nifi/N). The indices n and fi span the 
integers n = 0, . . . , A^ — 1 and fi — 1, . . . , (A' — 1) /2, respec- 
tively. 



(the last two relations are equivalent to the statement 
that Z commutes with all group elements, i.e., it is a 
central element). The set of relations ( p^ ) defines a cen- 
tral extension of the group Djq which we further denote 
ZDn. Extending Dn to ZDn by the central element Z 
enlarges the number of elements in the group by a factor 
of N , and the number of elements of ZDn is 2N'^ . 

In any particular irreducible representation, Z is rep- 
resented by a number (a A^-fold root of unity, since 
Z^ = 1): 

Z = exp(27r«Q/A^) (19) 

This formula establishes the relation between the repre- 
sentations of ZDn and the total charge Q on the array. 

Here we should mention that central extensions of 
symmetry groups are common in quantum mechanics. 
Indeed, the overall phase of the wave function has no 
physical meaning. Therefore the operator representing 
the product of two symmetry operations must equal the 
product of the two operators representing each of these 
operations only up to an overall phase factor. In other 
words, quantum mechanics admits not only linear rep- 
resentations of the symmetry grojip, but a more general 
class of projective representationsQ. At the same time any 
projective representation of a group corresponds to a lin- 
ear representation of a central extension of this group. 
Physical examples are numerous, including half-integer 
spin (projectively representing the rotation group), mag- 
netic translations (projectively representing geometric 
translations), and anyons (projectively representing the 
braid group) . An example resembling the analysis in-the 
present paper is given by Kalatsky and Pokrovskytil in 




(b) fc) 

FIG. 5. The irreducible representations of ZDn may be 
distinguished by the eigenvalues of the operator T. These 
eigenvalues must belong to the set of A-fold roots of unity 
{exp(27rm/A^)} (marked by filled dots on the unit circle). 
The dotted line is drawn to intersect the unit circle at points 
±\fZ. (a) TV odd: exactly one of the square roots ±\/Z is 
a A^-fold root of unity. It gives the eigenvalue of T involved 
in the one-dimensional representations D'* '" (marked by a 
small circle around the filled dot). The rest of allowed eigen- 
values combine in pairs to form two-dimensional representa- 
tions D*^"' (the corresponding pairs of dots are connected by 
dashed lines) . (b) and (c): the corresponding illustrations for 
even A and even Q, and for even A and odd Q respectively. 

their discussion of the spectrum of large spins in external 
crystal electric fields; the degeneracies then are described 
in terms of projective representations of a finite symme- 
try group which depend on the spin assuming integer or 
half-integer values. 

In the following we first review the irreducible repre- 
sentations of Dn and then discuss how it is modified 
when Dn is extended to ZDn- 

Consider first the case with an odd number of islands 
N . Using the defining relations (p^-(p^ we arrange the 
2N elements of the group Dn into the(iV -I- 3)/2 conju- 
gacy classes: {E}, {T^ ^-''}(^-^'/^ and {r™i?|^-i}. 
Firstly, we can construct two one-dimensional represen- 
tations D'-^'' where all operators are represented by num- 
bers. In these representations, the value of T is deter- 
mined uniquely: — 1, while R may take two values: 
D^' = ±1. In addition, we can find {N — l)/2 represen- 
tations I?'^' with dimensionality two. The explicit form 
of T and R in these representations is: 



6 



TABLE III. Character table of the dihedral group for A'^ even: the four one-dimensional representations are denoted by 
-D|j^'; the N/2 — 1 two-dimensional representations D*^' involve the roots = exp(27ri/i/A) with fj, = 1, . . . , N/2—1. 
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TABLE IV. Character table of ZDm for N even. The representations are classified by the charge Q which relates to the values 
of the central element Z — exp{2TviQ/N). The 2N one-dimensional representations D^ '" exist for even Q; the A'(A^ — l)/2 
two-dimensional representations denoted by D''' *" involve the roots = exp(27ri/i/A). The indices n, m, and fj, span the 
integers n = 0, . . . , A'^ — 1, m = 0, . . . , A''/2 — 1, and fi — 1, . . . , N/2 — 1, respectively; the index p. runs over the half-integers 
/i = l/2,3/2,...,(A-l)/2. 
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where the parameter ^ labeling representations takes in- 
teger values 1, . . . , {N— 1)/2, and = exp(27ri^/Af). We 
thus have found all the (A^-|-3)/2 irreducible representa- 
tions; they are listed in the character TABLE L 

The irreducible representations of the extended dihe- 
dral group ZDf4 are found in an analogous way. Assume 
N odd first. Then ZDn contains N{N + 3)/2 conju- 

gacy classes: {Z^^V, {Z^T- , Z-+-T-}'^Z~^;j:^i~^'''\ 
and {Z"-™T2™i? |^-J},^ro\ The construction of the 
irreducible representations again follows the scheme 
described above: for one-dimensional representations 
we find 27V solutions: D^*"' = \/Z and D*^*"' = 
±l/\fz, where we choose that branch of the square 
root which puts ^/Z onto one of the roots exp(27rm/A^) 
[for odd N either \/Z or — a/Z belongs to the set 
{exp(27rm/A^)|^Jo^}, see Fig. ||. The remaining two- 
dimensional representations may again be con- 
structed explicitly: 

= ^r' = (/i J). (21) 

where for ^/Z we again take one of the roots 
exp(27rin/iV), the parameter takes integer values 
1, . . . , (iV - l)/2, and = exp(27rj^/7V). With 2A^ one- 
dimensional and N{N — l)/2 two-dimensional represen- 
tations we then have constructed all irreducible repre- 
sentations of ZDn. The result is summarized in TABLE 
II. 

A similar analysis for the qubit loop with an even num- 
ber N of junctions produces the character TABLES III 
and IV for the groups and ZDn respectively. 

In summary, the energy levels of the Josephson junc- 
tion loop may be classified according to the irreducible 



representations of its symmetry group ZD^ at a given 
Q. The representation tables (Tables II and IV) look 
slightly different for odd and even N . All representa- 
tions are either one- or two-dimensional. Given an odd 
iV, for any value of Q, there are two one-dimensional 
representations and (A'^— l)/2 two-dimensional represen- 
tations. Given an even A^, for even Q we have four 
one-dimensional and {N/2 — 1) two-dimensional repre- 
sentations, while for odd Q there are no one-dimensional 
and A^/2 two-dimensional representations. The represen- 
tations may be distinguished by the eigenvalues of the 
operator T involved, see Fig. ^. 

Connecting back to our physical problem we can find 
the degeneracy of the states: for example, when A^ is 
even and Q is odd (see Fig. ^(c)) all states (including 
the ground state) are doubly degenerate. More work is 
needed to determine whether the ground state is degen- 
erate in the other cases. Inspection of the ground state 
wave functions in the limits Ej ^ Ec and Ej <C Ec 
shows that they transform with a representation involv- 
ing the eigenvalues T = 1 and T = exp{27TiQ / N) . In- 
deed, in the Ej ^ Ec limit, the state {(f>.i = —tt/N} 
is a T = 1 eigenstate, while the state represented by 
the points = tt/N, i ^ k; (j)k = tt/N — 27r}, 

k = 1, . . . , N hy virtue of ( p^ belongs to the eigenvalue 
T = exp(27riQ/Ar). In the limit Ej <t: E^ T measures 
the total momentum of the repulsive bosons which (for a 
diagonal capacitance matrix) equals the total momentum 
of the tight-binding fermions. The latter is easily read 
off Fig. ^ and accounting for the shift ttQ/N due to the 
boundary condition (originating from the half-quantum 
flux) we again obtain the two eigenvalues T = 1 and 
T = exp{27TiQ/N) for the two lowest states. Thus if 
Q (mod A^) the two lowest states combine into a 
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FIG. 6. The two lowest excitation energies for the sym- 
metric Josephson junction loop with A*' = 3 islands, diagonal 
capacitance matrix {C~^)ij = EcSij, and vanishing gate volt- 
ages Vi = 0. The ground-state energy at each value of Ec/Ej 
is subtracted, (a) in the Q = sector where the ground state 
is non-degenerate; (b) in the Q = 1 sector with the ground 
state doubly degenerate. Note the large (small) gap of order 
Ec (Ej/N) for the insulating state with Q = (the metallic 
state with Q = 1) in the limit Ec/Ej — » oo. 



two-dimensional representation (and the same in both 
hmits), resulting in a doubly degenerate ground state. 
For Q divisible by N , the ground state corresponds to 
a one-dimensional representation of ZDn, i.e., it is non- 
degenerate. As the ratio Ej /Ec is swept from one limit 
to the other, the ground state continuously evolves pre- 
serving its degeneracy, unless a level-crossing occurs with 
a level of different symmetry. We have confirmed numer- 
ically that such a level crossing does not occur for the 
cases N — ?, and iV = 4 using a diagonal capacitance ma- 
trix (C~^)ij, (the results for = 3 are shown in Fig. ||), 
and we believe that this property holds for any N > 3. 



VI. REQUIREMENTS FOR INTERFERENCE 
OBSERVATION 

The experimental observation of the quantum inter- 
ference effect discussed in this paper relies on several 
requirements on the Josephson junction loop. First, 
we have to assume that all quasi-particles are frozen 
out, which defines an upper bound on the temperature 
roughly estimated as T < T* ~ Age/ log(i^o^^sc), where 
1^0 is the density of states at the Fermi level, Age is the 
superconducting gap, and V is the volume of the sys- 
temE£l. For a micrometer-size aluminium loop similar to 
thai, used in the experimental work of van der Wal et 
ffl/.El, T* ~ Agc/15 ~ 100 mK and can be easily achieved. 

Second, we assume that no charge tunneling is possible 
onto the array. In practice, we may allow the array to 
be connected to an external reservoir via a large resis- 
tor (to be able to change the charge on the array by an 
overall shift of gate voltages). Its resistance then must 
be much larger than the characteristic resistance scale 
R* ~ (CA)~^, where C is the capacitance of the array, 
and A is the relevant energy scale (of the order of the 
splitting between the two lowest states). In the charge- 
dominated limit, A ^ Ec and R* is of the order of the 
resistance quantum h/Ae^, in agreement with the conven- 
tional conditioLLL-for charge quantization in the Coulomb 
blockade setupllEl. In the phase-dominated limit, A ^ Ec 
and R* is much larger than the resistance quantum. Also, 
under the condition Age > Eq°^ a 'weak' contact to the 
external world would allow an unpaired quasi-particle to 
escape from the loop, leaving only paired electrons in the 
system. 

Third, the islands and the junctions are assumed to 
be identical. The precision to which this symmetry has 
to hold in the phase-dominated limit may be simply es- 
timated from the condition that the tunneling actions 
agree up to small deviations of order one. The tunnel- 
ing action scales as S* {Ej/EcY^^ and therefore the 
required precision is {SEj/Ej), {6Ec/Ec) <C S~^. For a 
detectable splitting 5j-|must be not too large; for a qubit 
design, S is typicallyO'B taken to be of order 5-10, and 
these conditions may well be satisfied. 

The required precision on the magnitude of the exter- 
nal magnetic field (or, equivalently, on the flux through 
the loop) may be estimated from the condition that the 
level splitting due to the deviation of the flux from $0/2 
is much smaller than that from the tunneling between the 
qubit states. In the phase-dominated limit {Ej » Ec) 
this produces the condition i?j(S$/$o ^ ^ and hence 



In the opposite charge-dominated 



limit, it is sufficient to require that (5$/$o ^ 1- 

In conclusion we have discussed the spectral proper- 
ties of symmetric Josephson junction loops, devices sim- 
ilar to those recentl¥|-|proposed as potential qubits for 
quantum computingtrH. While in the charge-dominated 
regime with Ec ^ Ej the dependence of the ground- 
state degeneracy on the total charge on the island is a 
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simple charging effect, this degeneracy derives from a 
subtle quantum interference in the phase-dominated limit 
Ej ^ Ec ■ Several proposals on solid state realizations of 
quantum bits belong to the latter limit. From our anal- 
ysis it follows that the requirements for observing the 
charge dependence of the qubit-level splitting arc similar 
to those for the qubit operation (plus symmetry require- 
ments which are easy to satisfy). The interference effects 
studied in this paper then may serve as a good test of 
quantum coherence in such qubit designs. 

We thank M. Feigelman and O. Hallatschek for helpful 
discussions, and M. Troyer for providing the exact diago- 
nalization code. We thank the Swiss National Fonds for 
financial support. L. I. thanks ETH Ziirich for hospital- 
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